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The selection of a higher vapor deposition angle when growing a columnar thin film (CTF) leads to surface- 
wave propagation at a planar metal-CTF interface with phase velocity of lower magnitude and shorter 
propagation range. Acordingly, a higher angle of plane-wave incidence is required to excite that surface wave 
in a modified Kretschmann configuration. © Anita Publications. All rights reserved. 

1 Introduction 

An enormous body of literature on the propagation of electromagnetic waves localized to the 
planar interfaces of bulk metals and bulk dielectric materials has gathered during the last hundred 
years [1-3] . Such a wave is associated in quantum parlance with surface plasmon polaritons, resulting 
from the interaction of photons in the dielectric material and electrons in the metal. The quantum 
term is parsed as follows: the entities are localized to a surface; a plasmon is a collective excitation 
of electrons; and the dielectric material is polarized because of interaction with photons, thereby 
giving rise to the noun polariton. In the language of classical electromagnetics, the surface waves are 
p-polarized, not s~polarized. 

Although the dielectric material is usually considered to be isotropic and homogeneous, anisotropic 
dielectric materials (such as crystals) have been incorporated in surface-wave research as well [4,5]. 
As is known well, assemblies of parallel nanowires called columnar thin films (CTFs) can be used in 
optics in lieu of crystalline dielectric materials [6] . Provided the wave vector of an electromagnetic 
planar wave is oriented parallel to the morphologically significant plane of a CTF, a distinction be- 
tween p- and s-polarization states can be easily made [7]; hence, it can be conjectured immediately 
that surface waves can propagate on the planar interface of a homogeneous metal and a CTF. 

Surface-wave propagation (SWP) on a planar metal-CTF interface is bound to be affected by 
the morphology of the CTF. Columnar thin films fall under the general banner of biaxial dielectric 
materials for optical purposes. These thin films are grown by physical vapor deposition, whereby 
vapor from a source boat in an evacuated chamber is directed at angle Xv G (0,7r/2] towards a 
planar substrate, as shown in Figure [TJ Under the right conditions, parallel columns of the evaporant 
species grow on the substrate tilted at an angle Xv- The CTFs are composed of multimolecular 
clusters with ~ 3 nm diameter which, in turn, arc clustered in a fractal-like nature eventually 
forming columns with ^ 100-nm cross-sectional diameter, depending on the evaporant species and 
the deposition conditions. Once the evaporant species and the deposition conditions have been 
chosen, the vapor incidence angle Xv can be used to alter the morphology of a CTF significantly 
enough as to have optical consequences. Indeed, at visible frequencies and lower, a CTF may be 
regarded as a linear, orthorhombic continuum whose relative permittivity dyadic can be controlled 
by proper selection of Xv [8] • 

Our aim in this communication is to establish the influence of the CTF morphology, as captured 
by the vapor incidence angle, on SWP at a planar metal-CTF interface. Section 2 provides the 



derivation of the SWP wavenumber at a planar mctal-CTF interface, and Section 3 contains the 
solution of a boundary- value problem to excite a surface wave in a Kretschmann configuration [9,10] 
modified for practical considerations. An exp(— zwi) time-dependence is implicit, with uj denoting 
the angular frequency. The free-space wavenumber. the free-space wavelength, and the intrinsic 
impedance of free space are denoted by = uj^eg^a, \ = 27r/fco, and rja = \/ /io/co; respectively, 
with /io and Eq being the permeability and permittivity of free space. Vectors are in boldface, dyadics 
underlined twice; column vectors are in boldface and enclosed within square brackets, while matrixes 
are underlined twice and similarly bracketed. Cartesian unit vectors are identified as u^;. Uj, and u^. 
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Figure 1: Schematic of the growth of a columnar thin film. The vapor flux is directed at an angle Xd, 
whereas columns grow at an angle x ^ X"- The unit vectors Ut, Un, and ut, are also shown. 



2 SWP at Planar Metal CTF Interface 

Let the region z > be occupied by a CTF and the region z < by a metal of relative permittivity 
e„i. The relative permittivity dyadic of the CTF may be stated as [7] 

l^ty = EaUnUn +efcU^U^ +ecUfcUf,, (1) 

where ea,6,c are the principal relative permittivity scalars and the unit vectors 

u„ = -u^sinx + U;,cosx, = cosx + sinx , = -Uy . (2) 

All four constitutive quantities — ea,&,c and the column inclination angle x € (0,7r/2] — depend on 
the evaporant species and the vapor incidence angle Xv- The xz plane, containing the unit vectors 
u„ and u^, is the morphologically significant plane of the CTF. 

In order to preserve the independence of the p- and the s-polarization states, as stated in Sec- 
tion 1, the wave vectors in both half-spaces must not have a component orthogonal to the morpho- 
logically significant plane. Accordingly, the p-polarized electromagnetic field phasors in the metal 
half-space are given by 

E(r) = A„i [vi^ - l^u^^ exp {iko {kx - ig^z)] 
H(r) = A„i 77^7^ ^ Wy exp [ifco {xx ~ iq^z)] 

where Am is a complex-valued scalar of finite magnitude, qm = +V x'^ — ^m, and xkaUx may be 
regarded as the wave vector of the surface wave. We must have Re [qm] > for SWP. The p-polarized 



< 0, 



(3) 
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electromagnetic field phasors in the CTF half-space are given by [7] 



E(r) = A, 



H(r)=A,,;- 



exp [ika {xx + iqcz)] 



i^ge-(ea-e6)sinxcosx ^ 
>i2 — (e^ cos^ X+^fa sin-^ x) ^ 

-iqcjeg cos^ X+<^b sin^ x) + ^(<^a-<Si)) sin x cos x 
>t^ — (ca cos^ X+^b sin-^ x) 



Uj, exp [z/cq + iQcz)] 



2>0, (4) 



where Ac is a complex-valued scalar of finite magnitude, and Qc must be obtained by solving the 
quadratic equation 

>f^(eQ sin^ X + cos^ x) ~ "^-ixqci^a - £&) sin x cos x - qli^a cos^ x + £& sin^ x) ~ ^a^b . (5) 

We must choose Re [q^ > for localization of energy to the metal-CTF interface. Let us note 
parenthetically that does not enter our analysis. 

The boundary conditions at the interface z = lead to the two equations 



-iqcjig cos^ x+efc sin^ x) + ^(ea-eii) sin X cos X 
^ x'"^ — i^a cos'^ X+^b sin"^ x) 



4 H". 



which yield the dispersion equation [11] 



emX^ + iqmx{ea - e&)sinxcosx+ {qmqc - em)(eaCOS^x + ffesin^x) = 



(6) 



(7) 



for SWP. 



Equation cqn:SWP was solved numerically after choosing = —56 + 21i (typ., aluminum at 
An = 633 nm [2]) and choosing titanium oxide as the material of which the CTF is made. At 633- nm 
wavelength, empirical relationships have been determined for titanium-oxide CTFs by Hodgkinson 
et al. [8] as 

OH 2 



1.0443 + 2.7394 



1.6765 + 1.5649 



Xv 

7r/2 



Xv 

7r/2 



1.3697 



0.7825 



Xv 

tt/2 



Xv 

tt/2 



and 



tanx 



2.8818 tan x« , 



(8) 
(9) 

(10) 



where Xv and x are in radian. We must caution that the foregoing expressions are applicable to 
CTFs produced by one particular experimental apparatus, but may have to be modified for CTFs 
produced by other researchers on different apparatuses. 



Table 1. Dependence of >c on Xv for SWP at 633- nm wavelength 
on a planar interface of aluminum and a titanium-oxide CTF. 



Xv 




Im[>^] 


5° 


1.2619 


0.0179 


20° 


1.8506 


0.0174 


45° 


2.3244 


0.0326 


60° 


2.4688 


0.0410 


90° 


2.5783 


0.0488 
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Tabic 1 shows the computed values of the relative wavenumber >t for SWP on the chosen nietal- 
CTF interface. Clearly, as Xv increases towards 90°, both the real and the imaginary parts of yc 
increase. This means that an increase in the vapor incidence angle 

(i) reduces the magnitude of the phase velocity and 

(ii) increases the attenuation 

of the surface wave. Thus, a high value of Xv is inimical to long-range SWP. 



3 Modified Kretschmann Configuration 

In conformance with the Kretschmann configuration [9] for launching surface plasmon polaritons, 
the half-space z < is occupied by a homogeneous, isotropic, dielectric material described by the 
relative permittivity scalar t\. Dissipation in this material is considered to be negligible and its 
refractive index n\ = yjt\ is real-valued and positive. The laminar region < z < L,„ is occupied 
by a metal with relative permittivity scalar £,„ . A CTF of relative permittivity dyadic e occupies 

— Ct J 

the region Lm ^ z < — L,n + L^f- Finally, the half-space z > Ls is taken to be occupied by 
an isotropic, nondissipative dielectric material with relative permittivity €2 = n\> 0. This modified 
Kretschmann configuration is in accordance with practical considerations for launching surface waves 
[10]. 

A p-polarized plane wave, propagating in the half-space z < at an angle 9i G [0, tt) to the z axis 
in the xy plane, is incident on the metal-coated CTF. The electromagnetic field phasors associated 
with the incident plane wave are represented as [7] 



Ei,ic(r) = (— Uj- cos 6'i + U2 sin 6*1) exp + tiiz cos 6'i)] 
Hi„c(r) = —nirj^^ Uy exp [i{Kx + niz cos Oi)] 



z < 0, 



(11) 



where 

K = kaiii sin 9i . 
The reflected electromagnetic field phasors are expressed as 

Ere/(r) = TpivLx cos6'i + Uz sin 6*1 ) exp [i{Kx — nizcos6'i)] 



Hre/(r) = — fp '^1^0 ^ 6xp [i{nx — niz 



cos )] 



and the transmitted electromagnetic field phasors as 

Etr(r) = ip(-Uj, cos 6*2 + Uz sin6'2) exp{z [kx + 712 (z - ^s) cos 6*2]} 
Htr(r) = -tp n2rio^ v^y exp {i [kx + n2 (z - Ls) cos 6*2]} 

where 712 sin02 = ni smOi = x. 



z < 0, 



z> Ls 



(12) 



(13) 



(14) 



The reflection coefficient Vp and the transmission coefficient tp have to be determined by the 
solution of a boundary- value problem [7]. It suffices to state here that the matrix equation 



IK 



.-(2)1-1 



exp (i[£^j^]ict/ 



exp(j[P^jL„) .[^(1)] 



(15) 
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emerges, wherein the foUowing 2x2 matrixes are employed: 



(^)l _ 



- cos 9e 



cost 



^e{l,2} , 





K(ea-e6)smxcosx 
ea cos2 x+Cfa sin'' x 

-^f:n 9 ; : — 1 — 

cos- x+eb sm x 







'^/^o ujtoi^a cos2 x+ej, sin^ x) 
K(e(j — efa) sin x cos x 
ea cos^ X+^b sin-^ x 



(16) 
(17) 
(18) 



The solution of finaleq yields the reflection and transmission coefficients. The principle of conservation 
of energy mandates the constraint 



'n2 



, COSfc*! / 



(19) 



the inequality turning to an equality only in the absence of dissipation in the region < z < L-^. 
In order to study the excitation of a surface wave at the metal-CTF interface, the absorbance 



1 



/ Re [cos 02 
\ cos 6*1 



\tp\' 



(20) 



has to be plotted against the angle 9i. This was done for L,„ = 10 nm and Lctf = 1000 nm. As in 
Section 2, Cm = —56 + 21i was chosen, along with ea,b and x as specified by eqHl-eqH4. Whereas 
£1 = 9 was chosen to ensure that existence of a critical angle for total reflection in the absence of the 
10-nm-thick metal film, €2 = 5 was set for a lesser degree of constitutive contrast with the CTF. 

Figure 2 shows plots of Ap vs. 9i for the same five values of Xv as in Table 1. For comparison, 
plots of |rpp vs. 01 in the absence of the metal film arc also shown in this figure, in order to establish 
the critical (minimum) values of 9i for total reflection. For 9i greater than the critical angle for a 
specific Xv, we see very high Ap-peaks in Figure 2. 



Table 2. Dependence of 9i and x of the rightmost peaks in Figure 2 on x- 



Xv 


9i 


x ~ Til sin 01 


5° 


25.782° 


1.3037 


20° 


39.515° 


1.9088 


45° 


53.247° 


2.4037 


60° 


58.185° 


2.5493 


90° 


62.168° 


2.6530 



Values of 9i for the rightmost peaks of Ap against Xv in Figure 2 are provided in Table 2, along 
with the corresponding values of x ^ ni sin^i. The values of k in Table 2 are quite close to the 
real parts of those in Table 1, thereby indicating that the rightmost Ap-peaks in Figure 1 can be 
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Figure 2: Solid lines represent the absorbance Ap against the incidence incidence angle ^i, when ei = 9, 
—56 + 21i, Lm = 10 nm, Lctf = 1000 nm, and e^.b and x s^re specified by eqHl-eqH4. Dashed 
fines represent \rpf against Oi when Lm = 0. (a) Xv ~ 5°, (b) Xv ~ 20°, (c) Xv ~ 45°, (d) Xv = 60°, (e) 
X. =90°. 



attributed to the excitation of surface waves at the planar metal-CTF interface [2] . Not surprisingly, 
the value of 9i needed to excite a surface wave in the modified Kretschmann configuration increases 
as the vapor incidence angle Xv gets closer to 90°. 



4 Concluding Remarks 

Right from the time of Faraday [12], morphology is widely known to play a key role in the optical 
responses of thin films [13-15]. More recently, Hodgkinson and Wu made a comprehensive proposal 
to regard columnar thin films as equivalent to biaxial crystals [6], which functional analogy can be 
extended to sculptured thin films [7]. Therefore, it is appropriate that the excitation of surface waves 
at the planar interfaces of isotropic dielectric materials [10] and metals with CTFs be influenced by 
morphology. In this communication, we have shown that the selection of a higher vapor deposition 
angle when growing the CTF, with consequent influence on the morphology of the CTF, will lead to 
a surface wave at a planar mctal-CTF interface with phase velocity of lower magnitude and shorter 
propagation range; concomitantly, a higher incidence angle will be required to excite that surface 
wave in a modified Kretschmann configuration. 
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